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Abstract: We study the propagation of bubbles of new vacuum in a radially inhomogeneous 
Lemaitre-Tolman-Bondi background that includes a cosmological constant. This exemplifies 
the classical evolution of a tunneling bubble through a metastable state with curvature in- 
homogeneities, and will be relevant in the context of the Landscape. We demand that the 
matter profile in the LTB background satisfy the weak energy condition. For sample profiles 
that satisfy this restriction, we find that the evolution of the bubble (in terms of the physi- 
cally relevant coordinates intrinsic to the shell) is largely unaffected by the prsence of local 
inhomogeneities. Our setup should also be a useful toy model for capturing the effects of 
ambient inhomogeneities on an inflating region. 
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1. Introduction 

Although bubble propagation in homogeneous cosmological backgrounds has been amply 
studied [l]-[3], bubble propagation in inhomogeneous backgrounds is virtually unexplored 
[4]. Since tunneling is of enormous importance [5, 6] in populating the string Landscape 
[7, 8, 9, 10], we expect that the evolution of a new Universe through an old one, is worthy of 
study. 

One reason for the lack of emphasis on inhomogeneous backgrounds is the cosmic no-hair 
theorem [11, 12, 13]. This theorem states that the evolution of an expanding universe that 
contains a cosmological constant will eventually be dominated by vacuum energy, and will 
asymptote to a de Sitter space. It will imply in particular, that inhomogeneous backgrounds 
containing a cosmological constant (like the ones constructed on Landscape vacua) will evolve 
toward 'ironing out' such inhomogeneities, and therefore it will make sense to restrict the study 
of bubble propagation to homogeneous backgrounds. But fluctuations can still clump locally 
even in an asymptotically de Sitter space, and therefore the question of evolution of vacuum 
bubbles in an inhomogeneous spacetime is still relevant. In the context of eternal inflation 
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and the string landscape, we have to face up to the possibihty that the Universe that we 
inhabit is a result of bubble evolution through an inhomogeneous relic Universe. In principle 
this could give rise to experimental signatures in the form of CMB fluctuations, about which 
we comment in the final section. 

One particular work that stands out in this context is that of Wald [14], who did a 
comprehensive study of of the late time evolution of anisotropic but homogeneous spaces 
with a cosmological constant. In this restricted setting, he was able to define the conditions 
for the cosmic no hair theorem to be satisfied by taking advantage of the Bianchi classification 
of homogeneous, anisotropic metrics. Wald studied the time evolution of such metrics when 
the energy momentum tensor is the sum of two components: a cosmological constant and 
a term that satisfies the dominant and strong energy conditions. He found that with the 
exception of the Bianchi IX, the other models will always asymptote to a de Sitter space, 
within a time scale In the case of the Bianchi IX models, the future behavior depends 

on the relative sizes of the cosmological constant and the spatial-curvature. Only when the 
latter exceeds the cosmological constant will the time evolution not be asymptotically de 
Sitter. 

It is often claimed that inflation solves the horizon and flatness problems but, as already 
remarked in its early years, these claims arc based on the assumption that the pre- inflationary 
state of the universe can be described by a homogeneous and isotropic Fricdmann- Robertson- 
Walker (FRW) space. Under these strong assumptions, inflation does indeed reduce consider- 
ably the amount of flne tuning necessary to explain the current experimental measurements. 
In general, however, one expects a high degree of inhomogeneity in the initial conditions, with 
some regions of the universe expanding while others contract. Holographic cosmology might 
be the only exception to this rule because it naturally produces a homogeneous and almost 
flat space-time [15, 16, 17]. 

The question of whether an inflating region can continue to inflate if the ambient region 
is inhomogeneous, has been explored by several authors for different pre-inflationary back- 
grounds. In a set of papers [19, 18, 20], Goldwirth and Piran sought to answer the question 
of whether large inhomogeneity in the very early universe could prevent it from entering a 
period of inflation. In their work, inflation was not driven by a cosmological constant but by 
an inflaton fleld. Their analysis is numerical and restricted to an inhomogeneous but isotropic 
and closcd-imiverse background. They found that a large initial inhomogeneity could indeed 
suppress the onset of inflation, and that the inflaton fleld must have a suitable value over 
a region of several horizon sizes in order for inflation to start. This last flnding raises the 
issue of an acausal initial condition. Other numerical simulations disagree with these results 
claiming that inhomogeneous backgrounds will have enough inflation to explain current obser- 
vations [27, 28]. Analytical computations that use the long wave-length approximation seem 
to conflrm that large inhomogeneities of the spatial curvature prevent the onset of inflation 
[21, 22, 23]. 

In this work, we will not be able to conclusively settle this issue once and for all. Instead, 
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our goal here is to study a particular problem that we believe captures some of the same 
physics by playing local curvature effects against the might of the cosmological constant. 
The problem we consider is the classical evolution of a vacuum bubble in an inhomogeneous 
background, containing a cosmological constant and dust. This background, which we will 
refer to as the Lcmaitrc-Tolman-Bondi (LTB) space-time ^, will be the simplest possible 
inhomogeneous space, spherically symmetric and with only one center [24]-[26]. The bubble 
will be created in a region of space that is expanding but will encounter through its evolution 
regions of varying curvature. In a nutshell, the problem is that the critical size for bubble 
nucleation depends locally on the scale factor. As this factor changes as a function of the 
position of the shell, it is possible for a supercritical bubble to become subcritical ^. Our 
analysis will be confined to the evolution of the bubble, it will not address the interesting 
issue of tunneling probabilities in inhomogeneous backgrounds. This latter issue is very 
complicated, some attempts to address it have been made in [29]. 

Before embarking on describing the details, we emphasize that we are not dealing with the 
most general inhomogeneous bubble evolution imaginable. There are three major assumptions 
we make to render the problem tractable, and we list them below. The results of this paper 
should therefore be regarded as exploratory, rather than conclusive. 

• The first assumption we make is that the ambient spacctime into which the bubble 
evolves is the radially inhomogeneous LTB form. This allows us to work with an exact 
solution of Einstein's equation, but the restricted form of the metric and matter in this 
solution implies that we are not dealing with the most general radially inhomogeneous 
cosmology admissible. But as we explain in more detail in the next section, this still 
leaves us with a great deal of flexibility in choosing the inhomogeneity profiles. On the 
inside of course, we are dealing with a vacuum bubble, so we take the matter there to 
be only a cosmological constant without loss of generality. 

• The second assumption that we make is that the evolution of the new phase can be 
adequately captured using a thin-wall approximation. In particular this means that we 
can deal with all the transient phenomena between the two regions using a thin layer. 
The more general scenario would be to model the matter leakage from one region to 
the other using a gradual profile, but we believe the thin-wall captures the phenomena 
involved to a first approximation. 

• The above two assumptions together give us some traction. Once we make them, the 
problem that we need to solve reduces to that of solving the Israel junction conditions for 
a thin shell. But the junction conditions do not fully fix the evolution of the shell until we 
specify an equation of state on the bubble-wall. In the purely general relativistic context 

that we are working in, these equations of state are not fixed by the dynamics. So what 
we do here is to take the first step in this direction by considering a few equations of 

^Wc slightly generalize it from the original form to allow for a positive cosmological constant, 
bubble that can become subcritical in the future is a fluctuation, not a phase transition. 
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state of the form p = wp, for a range of values of w. We make some comments elsewhere 
on how one might try to fix this equation of state using field-theoretic arguments, but 
the full solution of this issue in the complicated background that we are working in, is 
likely to be difficult. In any event, as a start, we settle here for an exploration of these 
equations of state, leaving a more thorough analysis for future. 

We investigate a few different scenarios using our numerical formalism. The first is the ho- 
mogeneous limit of LTB, namely Robertson- Walker cosmologies with or without matter. We 
reproduce and extend known results in the literature using our methods, and offer some com- 
ments and new perspectives. Particular attention is paid to the case of bubble-propagation in 
an FRW universe with a Big Crunch in the future. We find that there are new ways (details 
are in a later section) in which the bubble can evolve in this case: it can even recoUapse. The 
final scenario wc consider is an LTB space that is truly inhomogeneous. In our examples, 
we pick our inhomogeneity profiles to have a fluctuation that is as sharp as possible while 
still respecting the weak energy condition. Surprisingly, we find that the qualitative bubble 
evolution (as seen in the physically relevant coordinates intrinsic to the shell or in the inside 
coordinates) is essentially the same with or without the fluctuation. This suggests that the 
energy conditions work towards the smoothcning out of exotic bubble evolution patterns. 

In the next section we shall describe the LTB model that represents the background on 
which the bubbles will propagate. In section 3, we will review the propagation of bubbles 
in a homogeneous background. We will consider the already known case where the outside 
background is simply a de Sitter space, either fiat or closed, and the novel case where matter 
is added to the cosmological constant. In this latter case, we shall explore the interesting 
situation of a bubble moving in a background that evolves from expanding to collapsing. 
In section 4 we shall study the truly inhomogeneous LTB background. Our results, future 
directions and cautionary remarks are put together in the Conclusion. We have relegated to 
the Appendix, an explanation of the known formalism to study the propagation of bubbles 
in general relativity. 

2. Setup: The Bubble and the Background 

A propagating bubble divides the space-time into three regions: outside, shell and inside. 
To give a proper description of this space-time in general relativity we will make use of the 
junction conditions, first presented by Israel [1] and extensively used since then by many 
authors. In particular, we will follow the implementation of this formalism developed by 
Berezin, Kuzmin and Tkachev [2]. 

The model we will study consists of a bubble of true vacuum propagating on a metastable 
state whose energy-momentum tensor contains a higher cosmological constant and dust. The 
bubble is assumed to be a thin-shell, with a perfect-fluid energy-momentum tensor whose 
equation of state will be allowed to vary. The true vacuum will be assumed homogeneous and 
isotropic while the outside background will be assumed to be spherically symmetric about 
one point in space. The line element of the outside region [30] is: 
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,2 ..2 [0'{t,r)+ra'{t,r)]'^ 2 2,. \ i-,r>, /oi\ 
as = at — — 2^ —dr — a [t, r)r dih (2.1) 

The expansion parameter a{t, r) is both a function of time and of the radial coordinate r. 
Partial derivatives with respect to the radial coordinate will be represented by a prime. The 
function i?(r) is an arbitrary function of the radial coordinate only and will be taken to be 
positive everywhere. We will study the evolution of the bubble for different choices of R{r). 
It can be seen from these equations that we can choose a function R{r) and then integrate 
for a{t, r) at each r. Then we can use this a{t, r) to define the d{t, r) in (2.4), with the caveat 
that the choice of R{r) has to be made in such a way that this energy density needs to be 
positive. Therefore, the description of the spacetime is essentially complete (apart from fixing 
certain initial conditions which we will get to in a minute), once we specify the function i?(r). 
Because it is our prejudice that curvature will interfere with the expansion of the bubble 
when its value is comparable to the cosmological constant, we will choose R{r) such that it 
can locally overwhelm the effect of the cosmological constant. 

The equations of motion for the expansion factor and the dust density d{t, r) , in units 
where SttG = c = 1, are: 



^ dta{t,r 
\ a{t,r) 



+ 



a^{t,r)R^{r) 



-d{t,r)a^{t,r){a{t,r) + ra! {t,r)) 



A 

a^{t, r) 
A 



A 



out 



(2.2) 
(2.3) 



^4 is a constant. These equations reduce to the familiar Friedmann-Robertson- Walker equa- 
tions in the limit where R{r) is independent of the radial coordinate. Besides choosing R{r) 
the solution to the equations of motion will involve the choice of the initial condition a(io, f). 
For simplicity, we will choose this function to be independent of the radial coordinate r. 

The results that will be presented in the forthcoming sections correspond to different 
choices of the function R{r). For different choices of R{r), the collapse- profile can be com- 
pletely different. One way to understand this is to notice that the Friedmann-type equation 
here can be thought of as an energy equation for a particle of zero total energy in a potential 

For those values of r for which the maximum of this potential is above zero, if the universe 
starts at a big-bang, it will recoUapse, but for other values of r there is no recoUapse. 

The bubble will be assumed to be a thin-shell (described by a hypersurface S) and 
spherically symmetric. Its energy momentum tensor will be assumed to be of a perfect-fluid 
type: Sr^ = a, Sg^ = S(j)'^ = —P, P = wa and with the metric: 

ds'^l^ = dr'^ - /9^(r)d02 (2.5) 
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The interior of the bubble will be assumed homogeneous and described by 



ds^ = dT^ -¥{T)[,^^-^ + z^(Kl2] (2.6) 



l + z 

with the following equation of motion for the scale factor: 

J^)^(t)»^(-) + ^ (2.7) 

We are assuming the inner space-time of the bubble to be open as derived from a tunneling 
process [31]. 

Given this setup, our aim is to understand the evolution of the bubble. The dynamics 
of the bubble is governed by the Israel junction conditions, the details of which are relegated 
to an appendix. The coupled differential equations that govern shell evolution turn out to be 
(7.8), (7.10), or in a more explicit form, (7.22), (7.24). 

As always in general relativity, along with the dynamical equations, we also need an 
equation of state to fully specify the dynamics. Determining the equation of state between 
a and P on the shell from first principles is a difficult problem that requires a field theoretic 
model for matter on either side. Since the LTB metric for a generic choice of the curvature 
profile R{r) is known only numerically, this is doomed from the start. In the absence of 
a detailed understanding of this phase transition we will confine our study to perfect-fluid 
shells, P = wa, and explore various values of w. 

It turns out that in the absence of dust on the outside, the shell-evolution equations 
can be written entirely in terms of p{t) because (7.10) becomes trivial, but this is no longer 
possible when dust is included. Indeed, equation (7.8) can be written as (reproduced from 
the appendix): 

= -1 + SV, (2.8) 

j2 -^in < ( ^ < ^ ( ^out ~ ^in A 



In general, to fix the evolution of B{t) we need (7.10), along with a knowledge of the scale 
factor through numerical solution of (2.2). 

We also get constraints on the parameters and the allowed phase space from the positivity 
of energy, in particular, the condition that a be positive. The explicit relations that we will 
use later can again be found in the appendix. 



3. Bubble Expansion in a Homogeneous Background 

The study of the bubble evolution can be done in different sets of coordinates. When dust 
is present, the problem of interest is most conveniently analyzed in the coordinates of the 
outside background. To familiarize ourselves with these coordinates we will devote this section 
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to analyze a simpler problem. The simplification will come from restricting the outside 
background to be homogenous (-R(r) = constant). First, we will redo the well studied motion 
of a vacuum shell in a cosmological constant background. Then we will add a dust energy 
density to the cosmological constant. 

3.1 Homogeneous Background without Matter 

The Israel junction conditions (7.8)-(7.10) tell us that it is consistent to have a vacuum shell 
(surface energy is constant and u; = — 1) in a transition that separates two de Sitter spaces 
with different values of the cosmological constant. Since this problem has been solved many 
times in the literature we will only quote the results [2] . 

In terms of the intrinsic bubble coordinates, the motion of the shell is given by: 

p{t) = — cosh Bt 
B 



where 



3 V4 3(7 

regardless of the curvature of the outside background. The effect of this curvature only be- 
comes manifest when we express the motion of the bubble in terms of the outside coordinates. 
In these coordinates the evolution of the position of the bubble is given by: 



dx _ -(1 - w)\Ji1 -Wm)^ - 1){B^ - f )(1 - g) 

dt a^B^ - ^) ^ • 

where the further assumption that Aj„ = and Aout = A has been made. The evolution of 
the scale factor is given by the well known expression 

da/dt\'^ 1 A 



a2i?2 3 

For times t » jj = -^3/A, (3.1) reveals the following asymptotic behavior : 



x{t) ^ i?sin _ 0{e-"') if i? > and finite (3.3) 



x{t) constant + 0(e"^*) if ^ oo (3.4) 




Figure 1: Time evolution of the bub- 
ble in the outside coordinates, x[t], for 
several initial sizes, Xinit = 20, 55, 90. 
The background is assumed homoge- 
neous (i? = 100) and without matter. 
The scales on the axes depend on the pa- 
rameter values we have chosen and are 
therefore not terribly important. 




Figure 2: Time evolution of the bub- 
ble in the outside coordinates, x[t], in ho- 
mogeneous backgrounds without matter. 
Colors correspond to different choices of 
R; i? = 100 red, R = 150 green, R = oo 
blue. 



The value of the constants in (3.3) and ( 3.4) is very dependent on initial conditions 
as can be seen in Figure 1. Also (3.5) corroborates the well known fact [31] that bubbles 
propagate asymptotically in time at the speed of light in the thin wall limit of A <C 3(T. 

A space with positive curvature only makes sense if a(0) > -j^. Once the initial value 
exceeds this critical value the space expands forever. The effect of the curvature becomes 
eventually negligible. Thus it is not surprising that the evolution of the bubble is similar 
regardless of the curvature. There is a difference in the possible asymptotic value that x can 
take; in the positive curvature case x always has to remain below R. The rate of the expansion 
is asymptotically given by B in the shell coordinates and by H in the outside coordinates. In 
terms of the inside coordinates the motion of the bubble is given by: 

z(r) = i^i + s2(r-ro)2 

This expression can be derived from the metric matching condition at the shell (7.3) and 
assuming that 6(To) = 1, which is a good assumption when Aj„ = 0. Tq is the inside time 
at which the bubble is created. For late times the bubble propagates at the speed of light as 
expected [31]. Ultimately, we are always interested in the motion of the bubble as viewed by 
the observer inside the bubble. The moral of this example is that the motion, as viewed by 
the inside observer, is the same regardless of the amount of curvature outside the bubble. 

In our equations above (and in the ones that follow) we have chosen the natural units 
c = 8itG = 1. Since we are dealing only with classical evolutions, h never comes up, so we have 
the freedom to pick one more unit. In the numerical evolution that we will undertake later, 
this means that there is rescaling of the coordinates (by a dimensionful parameter) that we 
are still free to do. We will choose the various parameters (like A, a, A) within a few orders of 
magnitude: the actual values themselves are not too relevant because of the above-mentioned 
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freedom in choosing the scale for the coordinates. PhysicaUy we expect the cosniological 
constant to reflect the scale corresponding to the minimum of some effective scalar potential, 
whereas the bubble stress-tensor should reflect the potential barrier between the two minima. 
For a generic potential, we expect that the scales involved will be not too different in order 
of magnitudes. As far as the evolution of the shell goes, the only place where we could have 
trans-Planckian effects becoming significant is when there is a singularity in the spacetime 
itself. In particular the fact that the shell surface is an actual physical wall, and not a surface 
at conformal infinity (like a horizon) , means that we do not have to worry about finite energy 
modes observed far away being red-shifted trans-Planckian modes [32, 33, 34, 35]. 

The effects of curvature on the bubble propagation, as seen by the outside observer, are 
illustrated in Figures 1, 2 and 3. To simulate these evolutions we have chosen the following 
values for the parameters, in units of SttG = 1: 

A = 3 X 10~^ 

If out ~ ^in = +1 

(Tinit = 10-3 (3.6) 

These values for the parameters can not be chosen independently since they have to satisfy 
the condition (7.28) that yields: 

a < 2^ = 6.3x10-^ (3.7) 
Also from the equation (3.1) the lower and the upper bound on x are: 
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1 

oB 



< x< R 



(3.8) 



We have chosen the initial size of the bubble in accordance with this range. 

In particular, Figure 1 shows the evolution of the bubble in the outside coordinates 
for various initial sizes. Note that the maximum size of the bubble in these coordinates 
corresponds to the limit x = R. 

Figure 2 shows the effects of curvature on bubble propagation by comparing evolutions 
of bubbles with the same initial size Xinit, but in different curvature background. We see that 
more curvature corresponds to a slower evolution of the bubble in the outside coordinates; 
however, the corresponding bubble evolution p[r] in the bubble coordinates is exactly the 
same, irrespective of the amount of background curvature, and it is shown in Figure 3. 



o 

0. 001 




Figure 4: Time evolution of the surface energy density a[t\ on the bubble. The background is 
assumed homogeneous {R = 100) and without matter. Colors correspond to different equations of 
state; w = 1/3 red, w — Q yellow, w = —1/3 green, w = — 1 blue. 



In addition, one might wonder what would happen if we did not assume the bubble to 
be a vacuum bubble, namely \l w ^ —1. Then the evolution of the energy density a on 
the bubble is governed by (7.24). In the absence of dust on the outside the solution to this 
equation is: 



c 



p2(l+w) 

where ^ is a constant. When substituting this solution on (2. 
equation for pij) 



(3.9) 

we obtain the following 



-1 + 



c 



+ 



A 



4p(i+2«,) 3^ 



P 



3+2to 



(3.10) 
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Figure 5: Time evolution of the bub- 
ble in the coordinates of the bubble, 
p[t]. The background is assumed homo- 
geneous (R — 100) and without matter. 
Colors correspond to different equations 
of state; w = 1/3 red, w = yellow, 
w = —1/3 green, w ~ —1 blue. 



Figure 6: Time evolution of the bub- 
ble in the outside coordinates, x[t]. 
The background is assumed homoge- 
neous {R = 100) and without matter. 
Xinit = 20. Colors correspond to dif- 
ferent equations of state; w = 1/3 red, 
w ~ yellow, w = —1/3 green, w = — 1 
blue. 



For p large, the solution to this equation takes the form: 



P{r) ~ - 



(3.11) 



l-(l + 27x;)Ap(ro)i+2-(r-To) 



l/(l+2«)) 



This analysis of the asymptotic behavior is also captured in the results plotted in Figure 4, 
and Figure 5. 

Figure 6 depicts the corresponding bubble evolution in the outside coordinates, x[t]. 

In summary, in this section we have considered a closed FRW background that will 
eventually asymptote to a de Sitter space. The outside curvature doesn't affect the evolution 
of the vacuum bubble when it is described in terms of the shell coordinates, and therefore, 
nor in terms of the inside coordinates. The curvature only makes a difference when the 
evolution is studied in terms of the outside coordinates. In these coordinates, the effect of 
a bigger curvature is to slow down the evolution of the comoving bubble coordinate. The 
asymptotic value of this coordinate depends on the initial condition but is always smaller 
that the curvature radius. When the equation of state of the shell is not that of vacuum, the 
evolution in the shell coordinates is much faster than that of vacuum but is also independent 
of the outside curvature. When described from the point of view of the outside coordinates, 
this evolution is faster than the corresponding one for the vacuum shell, but is qualitatively 
similar. 

3.2 Homogeneous Background with Matter 

To study the evolution of the bubble in a space that might contract we have to modify the 
outside energy density. One option is to introduce matter. In this instance the evolution of 
the scale parameter is given by 
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(ia\2 1 A 2 ^ ^ 

dt) -a=^ 

A < ^ and a(0) is small enough the expansion of the universe will reverse into 
contraction. With this new form of outside energy density, however, it is no longer consistent 
with the matching conditions to have a constant surface energy density on the bubble, not 
even when a + P = 0. The evolution will be given by the equations: 



(3.12) 



_ - (l - ^) ^ ± V (l - f ) - 1) - (t + ^)) 

di ~ xa2(S2 - ^) 

$ = -2 f + '-^) (a + P) + 7o.. , —'4 (3.13) 

V « ^ J ^-(dxM)2a2 + l-f 

a 1 /A yl 
4 V 3 ^ 



5=T + -(o+-3) (3-14) 



As in the previous section in order to solve these equations an assumption about the equation 
of state on the bubble was needed. We did make the assumption that the shell is made of a 
perfect fluid with equation of state P = wa, and explored several values for w. For the dust 
energy to be comparable to the cosmological constant at early times, we choose 

A = IQ-^ (3.15) 

in Planck units. The critical curvature which is the minimal curvature needed for the space 
to turn around and eventually collapse is given by: 

/9AA2\ "^/^ 

Rcr = [-^) = (3.16) 
The position of the maximum of the potential V[a] is given by: 

/3^\ 

amax = i^-^j = 1.7 (3.17) 

In our simulations we choose the initial value of the scale factor a to be Uinu = a(0) = 1. 

The outside background geometry will be completely specified once we choose the amount 
of curvature. We will consider several cases. 

3.2.1 R > Rcr 

The universe will always expand. If at the instant the bubble is created a(0) is to the left of 
max of the potential (2.4), then the universe will experience a period of slower growth until 
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it eventually goes over the max and the expansion becomes dominated by the cosmological 
constant. 

We will now investigate bubble propagation on such background. We are again assuming 
that the inside region has zero energy density, and that ^cmt = 7m = +1- In this case, the 
condition (7.28) yields: 



.<2^| + ^ (3.18) 

This constraint will be the strongest as a — oo, when, as is the case without matter, it 
reduces to: 

(7 < 6.3x10-2 (3.19) 

Once again, we choose: 

(Tinu = 10-2 (3.20) 

where ainit is given in Planck units. Prom the equation (3.12) the lower and the upper bound 
on X are: 

— <x<R (3.21) 
aB - - ^ ' 

Wc will choose the initial size of the bubble in accordance with this range and compare 
the evolution for several equations of state on the surface of the bubble. Decreasing R (and 
therefore increasing the curvature of the space) will have the effect of slowing down the 
evolution, as long as i? > Rcr- 

Pirst let us consider the case with very little curvature, R = 500. Pigure 7 shows the time 
evolution of bubbles of several initial sizes in the outside coordinates. In most cases bubbles 
grow without reaching the upper bound, x = R. Smaller w slows down the expansion rate, 
and in particular, the choice of w; = — 1 soon leads to a contracting bubble in the outside 
coordinates. Depending on the initial size of the bubble, w = —1 bubbles either reach the 
lower bound, x = which results in the breakdown of the simulation, or eventually stabilize 
and stop contracting. 

Time evolution in the coordinates on the bubble (Figure 8) reveals behavior similar to 
the case without matter. Namely, for the bubbles that do survive, the evolution is slower 
than in the case without matter, but the overall qualitative behavior stays the same. 

Next we turn to the case with more curvature, but still such that R > Rcr- Prom Pigure 
9 we see that in most cases bubbles again grow, but now they asymptote to the upper bound, 
x = R. 

The corresponding bubble evolution in the coordinates on the bubble is shown in Pigure 
10. The overall effect of more curvature is a further slowdown in the time evolution of the 
bubbles. 
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Figure 7: Time evolution of the bub- 
ble in the outside coordinates, x[t], for 
several initial sizes, Xinu = 10,50,100. 
The background contains homogeneous 
curvature {R = 500) and matter {A = 
10"''). Colors correspond to different 
equations of state; w = 1/3 red, w = 
yellow, w = ^1/3 green, w = —1 blue. 



Figure 8: Time evolution of the bubble 
in the bubble coordinates, p[t], for the 
same background conditions as in Figure 
7. Colors correspond to different equa- 
tions of state; w — 1/3 red, w — yellow, 
w = — 1/3 green, w — ~l blue. 



Figure 9: Time evolution of the bub- 
ble in the outside coordinates, x[t\, for 
several initial sizes, Xmit = 10,50,100. 
The background contains homogeneous 
curvature (i? = 120) and matter {A — 
10"''). Colors correspond to different 
equations of state; w = 1/3 red, w = 
yellow, w = —1/3 green, w = —1 blue. 



Figure 10: Time evolution of the bub- 
ble in the bubble coordinates, p[t], for the 
same background conditions as in Figure 
9. Colors correspond to different equa- 
tions of state; w = 1/3 red, w — yellow, 
w = — 1/3 green, w = —1 blue. 



In summary, in this section we have considered a closed FRW background that contains 
dust as well as a cosmological constant but that will eventually asymptote to a de Sitter 
space. In this case it is no longer consistent to have vacuum shells separating de interior and 
exterior region ^. We have studied the motion for different values of w. In shell coordinates, 

^The point of this section is to consider bubble creation for times sufficiently small for the dust density to 
make a difference. For long times the dust density will redshift and this background will become indistinguish- 
able from the previous one. 
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for w = —1, we have found that some bubbles do eventually become subcritical or stabilize 
and stop contracting. 



3.2.2 R<R, 
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Figure 11: Time evolution of the scale factor a[t] for R = 100. (A = 3 x 10'^, A = 10"^.) 

For curvature greater than the critical value (i.e. R < Rcr) the space will eventually 
collapse. To investigate this case we choose R = 100. The evolution of the scale factor is 
shown in Figure 11. 

In the outside coordinates, the time evolution of most bubbles will once again lead to the 
upper bound, x = R (Figure 12). Bubbles with w = —1 contract and hit the lower bound 
^ ~ a5' which stops their further evolution. 

Figure 13 reveals corresponding behavior in the coordinates on the bubble. Bubbles 
which do not hit the bound x = eventually collapse along with the collapse of the space 



In summary, in this section we have considered a closed FRW background that contains 
dust and a cosmological constant and that will eventually contract. As in the previous case it 
is not longer consistent to have vacuum shells separating the inside and outside region. De- 
pending on the initial size, bubbles with smaller value of w become subcritical and disappear 
while the others continue to grow and eventually collapse along with the space itself. 

4. Bubble Expansion in Inhomogeneous Backgrounds 
4.1 Generating a Curvature Profile 

In order to generate an inhomogeneous curvature profile, we choose: 



itself. 




(4.1) 
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X 




Figure 12: Time evolution of the 
bubble in the outside coordinates, x[t], 
for crunching background (Figure 11). 
Xinit — 10. R ~ 100. Colors correspond 
to different equations of state: w — 1/3 
red, 1/7 = blue. 



Figure 13: Time evolution of the 
bubble in the bubble coordinates, p[t], 
for crunching background (Figure 11). 
Xinit — 10. R — 100. Colors correspond 
to different equations of state; w — 1/3 
red, w = blue. 



in the definition of the LTB metric. Choosing (3 > 1 will allow us to generate a sharp drop in 
the function R{r), for the cases where the minus sign is chosen. The position and the width 
of the extremum will be regulated by the parameter 6, and the depth/height by 7. 

Furthermore, we want to make sure to satisfy the weak energy condition, namely that 
the matter density stays positive: 



This means that a(r) should nowhere fall faster than 1/r. 

It should be noted that (4.1) incorporates an approximation to a delta function in the 
curvature profile. Since a constant value of R{r) can be re-absorbed by a variable redefinition 
in the metric, we are interested in the physics associated with a spatially varying R{r). In 
regions of spacetime where R{r) is roughly constant, we expect that the evolution of the 
bubble would be as in the FRW case^. In the next sub-section, two representative examples 
which capture the essential features of these gradient effects are presented. 

4.2 Examples 

As an example, we first choose the minus sign in R{r), with a = 0.1, (3 = 3 and S = 20. In 
order for R to stay above Rcr, in this case we need 7 > 1.5. In addition, numerical simulations 
show that to satisfy the weak energy condition, 7 > 2.9, and since we'd like to investigate the 
sharpest possible profile, we choose 7 = 2.9. 

*A crucial ingredient in coming to this conclusion is the fact that the LTB scale factor can be solved as an 
ordinary differential equation in time t, point by point in the spatial grid. See Section 2. 




)) 



(4.2) 



a{t, r) + ra{t, r) > 
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Figure 14: The position dependent 
scale factor a{r, t) of the the LTB space- 
time, for the curvature profile depicted in 
red in Figure 15. 



Figure 15: R{r) works as an inhonio- 
geneity profile. Two choices are shown 
here, the red curve results in an inhomo- 
geneous spacetime. 




Figure 16: Evolution of bubbles in out- 
side coordinates, x{t), for the two curva- 
ture profiles shown in Figure 15. 



Figure 17: Evolution of bubbles in bub- 
ble coordinates, p{t), for the two curva- 
ture profiles shown in Figure 15. 



This choice of the function R{r) generates a profile in the evolution of the scale factor 
a(r, t), as shown in Figure 14. We will compare evolution of a bubble in such background with 
evolution in a curvature background without a fluctuation in the profile. The two choices of 
R(r) are shown in Figure 15. The resulting bubble evolutions are shown in Figures 16 and 
17. We see that even though in the outside coordinates bubble evolution is greatly affected 
by the curvature profile, such effect is absent for the bubble evolution in the coordinates on 
the bubble. 

For comparison, we will also choose i?(r) with the plus sign, and a = 0.1, /? = 2,7 = 0.7 
and 5 = 20, depicted in Figure 19. This results in a profile for a(r, t) shown in Figure 18. 
The corresponding bubble evolution is shown in Figures 20 and 21. Again, we find that the 
fluctuation does not qualitatively alter the bubble evolution. 

5. Conclusions 

We have examined the evolution of bubbles of true vacuum in several backgrounds (that 
always include a cosmological constant) to assess the effect of curvature on the propagation. 
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Figure 18: a{r, t) for curvature profile 
depicted in red in Figure 19 



Figure 19: Second example of an inho- 
mogeneous background: Two choices of 
R{r) are shown. 




Figure 20: Evolution of bubbles in the 
outside coordinates, x{t), for the two cur- 
vature profiles shown in Figure 19. 



Figure 21: Evolution of bubbles in the 
shell coordinates, p(t), for the two cur- 
vature profiles shown in Figure 19. 



The bubble evolution changes qualitatively when looked at in terms of the outside coordinates, 
but when studied in terms of the more physically relevant coordinates on the shell or its 
interior, the local curvature is largely irrelevant. In the presence of matter, we observe that 
some bubbles, with an equation of state close to w = —1, disappear. Those that survive do 
behave qualitatively the same. 

The explanation for this insensitivity to inhomogeneities is to be found in the weak 
energy condition: our background spacetime is consistent only when (4.2) is positive, which 
restricts the sharpness of the inhomogeneity profiles that are allowed. This could of course 
be a restriction arising from the ansatz for the background that we are working with. But it 
should be mentioned that in the work of Wald [14] that was mentioned in the introduction, 
the energy conditions played a crucial role in the proof of the cosmic no hair theorem for 
anisotropic (but homogeneous) cosmologies. 

We view this work as an attempt to test the requirements to arrive at a universe like 
ours without making the simplifying assumptions of starting with an FRW metric. There 
are many directions in which the work here can be extended. We have examined only a 
few inhomogeneity profiles here, a more thorough understanding of the evolution curves in 
more generic situations is certainly of interest. For all the various initial conditions and 
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collapse/inhomogeneity profiles that are allowed, the question of bubble evolution should 
be tractable with our machinery^. If we think of our setup as a toy model for inflation, 
it will definitely be of interest to study the genericity of inflation in the space of possible 
inhomogeneity profiles. Another important direction to follow is to evaluate the tunneling 
amplitudes for inhomogeneous ambient spaces. This is a notoriously difficult problem on 
which any enlightenment will be welcome. 

Another problem that is of interest in this context is whether there are signatures left on 
the new (inside) Universe due to the propagation of the bubble through the inhomogeneous 
ambience. The bubble can be thought of as a Casimir cavity with a moving wall and the 
cavity radiation (the analogue of the CMB in this context) will carry imprints of the external 
inhomogeneities. Because of the spherical symmetry, this is effectively a two-dimensional 
QFT problem, with moving mirror boundary conditions [40]. 
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7. Appendix 

In this appendix, we will derive the coupled differential equations that form the starting point 
for our numerical shell evolution plots. These equations are a direct consequence of Israel's 
junctions conditions^, which relate the discontinuity in the extrinsic curvature across the 
shell surface to the energy-momentum layer on the shell. The general expressions for these 
equations for the case of spherically symmetric junctions can be found in, e.g., [2]. Here we 
discuss a few of the ingredients that go into our specific problem by doing some illustrative 
computations ab initio. The aim is to stress some issues like choices of signs and coordinates 
which are often not adequately addressed in the litearture. 

We will work with the mostly minus {-|-, — , — , — } metric, and the signs on the energy- 
momentum tensor for the ideal fluid are fixed by 

= {€+p)UaU^ -pSj'. (7.1) 

^Though such a scan would admittedly be much more numerically challenging than the few cases we have 
undertaken here. 

® A pedagogical derivation of the junction conditions based on a distributional approach to Einstein's equa- 
tions can be found in [36]. 
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Einstein's equation takes the form 0^/3 = T^p, in natural units (SvrG = c = 1). Later, when 
we work with quantities defined on the shell, we will have analogous sign conventions for 
them as well. Throughout, Latin indices denote 3-dimensional objects defined on the shell 
hypersurface, Greek indices stand for 4-dimensional quantities, and semi-colon is shorthand 
for the covariant derivative. We will also need e" = ^j-, which are projectors (pull-backs) 
that can be used to project a 4-quantity on to the 3-surface. 

The first of Israel's two conditions says that the metric induced on the shell from the 
bulk 4^metrics on either side should match, and be equal to the 3-metric on the shell. The 
assumption of spherical symmetry restricts the form of the intrinsic metric on the shell to the 
form, 

ds3^ = dr^ - p{Tfdn2, (7.2) 

where r is the only independent coordinate on the shell, which we take to be the shell proper 
time. So looked at from the outside, on the shell, we can parameterize the coordinates as 
r = x{t) and t = t^r), and from the inside, T = T{t) and z = z{t). Since the metric from 
either side on the shell should agree with the 3-metric on the shell, we get two conditions 
from the inside (open-FRW), 

and two from the outside (LTB), 

a{t,x)x = p{r), [-) =l+ i_,2/^(,)2 UJ • (7-4) 

All the variables in these equations are thought of as functions of r. 

Now we turn to the second junction conditions, which determine the dynamics of the 
shell. These are expressed in terms of the extrinsic curvature: 

Kab = ria-p ef . (7.5) 

Here is the outward normal to the surface under question (since we are working with a 
closed shell, there is no ambiguity in making this choice). There are many different choices 
and sign conventions for the extrinsic curvature in the literature, we have made the above 
definitions so that the extrinsic curvature of a 2-sphere is positive. In particular, this differs 
from [2], but is in agreement with most of the other references/reviews on the subject [36, 
37, 38, 39]. With these prescriptions, the second junction condition takes the form 

[Kah\ - hab[K\ = Sab, (7.6) 

with hab denoting the shell metric, which in our case is (7.2). The 3-tensor Sab stands for the 
surface energy-momentum tensor, and we will assume it to have a perfect-fluid form analogous 
to (7.1). Square brackets stand for discontinuities across the shell: [K] = Kgut — I^in- 
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It can be shown [2], that the junctions conditions imply a conservation law of the form 

Sa';b+[^aTinp]=0, (7.7) 

where T*^/ , s is the bulk enere-y-momentum tensor on either side. The advantage of this 

a(out/m) ^ 

equation is that it is first order, and one can exchange^ one of the second order equations 
arising from the junction conditions with (7.7). This is indeed what we will do. 

Because of spherical symmetry and the form of the three metric, the extrinsic curvature 
KJ' has independent components K^'^ and Kg^ = K^, while the surface energy tensor 
contains S/ = a and Sq^ = = —P. So the independent relations that arise from the 
second junction conditions are 

- f = [Kl], (7.8) 
P=[K;] + [kI]. (7.9) 

Both a and P are purely functions of r by spherical symmetry. The conservation law (7.7) 
takes the form, 

^ + -^(a + P) + [T-] = 0, (7.10) 

where [T"] = [e'^Tanp]. The evolution of the shell is completely determined by (7.8) and 
(7.10), so our task then is to write down these differential equations exphcitly so that we can 
proceed with the numerics. 

We have different sets of coordinates in each of the three regions, but not all of these 
coordinates can be explicitly written in terms of the others. So we will write the matching 
conditions in terms of the LTB coordinates. The LTB metric is known only numerically (see 
(2.2), (2.4)), so there is no hope of writing it in terms of the other coordinates, thereby making 
this choice inevitable. Since the evolution is best understood in the shell coordinates, once 
we have the evolution curves in LTB coordinates, we will translate them numerically to p{t). 

We start with (7.8). To calculate the extrinsic curvatures, we need the normal vectors 
in the corresponding coordinates. We will start with the LTB side where the coordinates are 
x° = {t,x,9,(f)). The projectors are: 

^ / dt dx \ 

e^ = (0,0,l,0), e^ = (0,0,0,1). 

Since is the 4- velocity of the bubble, the normal np will be determined (upto a sign) by the 
two conditions u'^ria = and n^ria = — 1, where raising and lowering are done with the LTB 
metric. The second condition arises because our shell is timelike: a timelike shell is defined 

'^This is analogous to the fact that the covariant energy conservation equation is an integrabiUty condition 
for Einstein's equations. 
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by a spacelike normal. A simple calculation yields, 

(—(ax),xX (ax),xi \ 
, " =, / " ,0,0 (7.11) 
v/(l - x2/i?2) V(l - a^V-R') / 

where the dots are with respect to r. The choice of sign for the normal is encoded in jout, 
and is fixed by the condition that the normal should point from the inside to the outside. For 
an expanding shell, this means that 'jout = +1- For a collapsing shell 'jout = — 1- 
Using these, the Kg on the LTB side turns out to be, 

T^e _uee^ - fr, r" r, - ^ / (aa:),^ (ax) ± + jl ~ R'^)i \ 

Keiout) - h ^e;e - ^ [ne,e - r^^na j - -iou\ ^^^^ _ ^2/^2) ) " ^^■^^> 

We have used (7.4) to do sonic of the simplifications, and (ax),^ stands for {xdta) because 
at this stage x and t are unrelated variables: when we eliminate r to write x as a function 
of t, this will no longer be the case. Repeating the above calculation for Kq on the inside 
(open-FRW), with coordinates x° = {T,z,6,4>), we have 

Again, 7j„ = +1 when the radius of the inside region is increasing; otherwise, = —1. The 
junction condition becomes, 

'{ax),tiax),xX+{l-x'^/R'^)i\ / zh§^z + {I + z^)f\ _ a 



'Yout I 



By using (7.4), it is possible to rewrite this equation in a more standard form as: 

where Agut = -(1 - x^/R^) + {ax),t^, and Aj„ = -(1 + z^) + z'^{dh/dTf. The quickest way 
to demonstrate this is to start with the final expressions. For instance, for the LTB piece, we 
can expand p as {ax)^tt + {ax), the relation 

once, assemble a perfect square from the pieces, and then use (7.16) again to end up with 
the first piece in (7.14). Equation (7.16) here is just a rewriting of (7.4). An analogous 
transformation can be done for the open-FRW part of the equation. 

The advantage of the form (7.15) is that now we can use the LTB equation of motion 
(resp. the open-FRW equation of motion) to simplify the two terms further to write 

Ao«* = -l+(^ + ^)p^ (7.17) 
Ain = -l + ^p'. (7.18) 
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Now, we take up the task of writing these equations purely in the LTB coordinates, for 
otherwise, despite being correct, they will be of no practical value in computations because 
of the mixing of coordinates. Using (7.17) and (7.18), we can write (7.15) as 

f? = p^B^ - 1, (7.19) 

where \ \ a ~ 

t> = — — V 



i + 3 +y)) • (7-20) 



3 

This form is useful because we just have to focus on p and p because everything else is already 
manifestly in LTB coordinates. The idea is that instead of choosing p(r) as the curve for the 
bubble-evolution, we want to parameterize it as x(t). Since x and t are dependent variables 
on the shell, this is legitimate. So we write p = ax, and rewrite p as 



[xdta^ {ax),^ 



(7.21) 



Using this in (7.19), and using the LTB equation of motion for dta, we end up with. 



dx 



±x{(l - ^){a'x^B^ - 1) (a2i?2 - f - ^) } 



1/2 



(7.22) 



dt {ax),^{a?x'^B'^ -x'^/B?) 

This is the first of the two shell evolution equations (the other being (7.10)) in a form that 
is directly applicable for numerical simulations. The sign on the square root comes from the 
square root of (7.19) and has to be chosen according to whether the shell is expanding or 
contracting. 

The second of the coupled shell-equations is easily translated into the LTB coordinates 
as well. On the LTB side, matter takes the form dust -|- cosmological constant (2.2). So, 
using (7.11), 



T" -r^T*. I r-T-. - 7o«K«x),.fd(x,t)Vl-xVi^2 

^r(ont)-Mt«t + e.i.n.- i_a:yR^-{ax)l{dx/dt)^ ' "^^ 



On the inside, the only matter is the vacuum energy, which gives rise to T^^.^-^ = 0. These 
along with (7.4) can be used to write (7.10) in the form 

^ ^ (x^a+M^^^ ^ _ 7o..(ax),.td(x,t) ^ ^ 

^1 - f - {ax)%{dx/dtr 

Together with an equation of state relating a to P, and the LTB equation of motion (2.2), 
(7.22, 7.24) comprise the starting point for the various special cases we study in the main 
text of this paper. 
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When solving for shell dynamics, we should demand that the positive energy condition 
be satisfied, that is 

(7 > 0. (7.25) 
Using the formulas derived here (in particular (7.15)), we can write : 

^out - ^in = ^ + loutP<y{p^ - Ao„t)l/^ (7.26) 

which together with (7.25) implies: 

Ao«t - Aj„ > — — if 7o„t = +1 



Furthermore, defining: 



Aout — Aj„ < — — if ^out = —1 



and using (7.15), the following relationship holds: 

7m I e+1 I -7o«d 1=2 (7.28) 

Depending on the outer and inner geometry, and given the energy densities both inside 
and outside the bubble, there will only be a certain range of possible values for the energy 
density a on the surface of the bubble consistent with (7.28). We will use this to constrain 
the value of a in the various cases addressed in this paper. 
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